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ABSTRACT 

Massive black holes (MBHs) in contrast to stellar mass black holes are expected to substantially 
change their properties over their lifetime. MBH masses increase by several order of magnitude over 
the Hubble time, as illustrated by Soltan's argument. MBH spins also must evolve through the series of 
accretion and mergers events that grow the MBH's masses. We present a simple model that traces the 
joint evolution of MBH masses and spins across cosmic time. Our model includes MBH-MBH mergers, 
merger-driven gas accretion, stochastic fueling of MBHs through molecular cloud capture, and a basic 
implementation of accretion of recycled gas. This approach aims at improving the modeling of low- 
redshift MBHs and AGN, whose properties can be more easily estimated observationally. Despite 
the simplicity of the model, it captures well the global evolution of the MBH population from z ~ 6 
to today. Under our assumptions, we find that the typical spin and radiative efficiency of MBHs 
decrease with cosmic time because of the higher incidence of stochastic processes in gas-rich galaxies 
and MBH-MBH mergers in gas-poor galaxies. At z = the spin distribution in gas-poor galaxies 
peaks at spins 0.4 — 0.8, and it is not strongly mass dependent. MBHs in gas-rich galaxies have a 
more complex evolution, with low-mass MBHs at low redshift having low spins, and spins increasing 
at larger masses and redshifts. We also find that at z > 1 MBH spins are on average highest in high 
luminosity AGN, while at lower redshifts these differences disappear. 
Subject headings: Black hole physics — Galaxies: active, nuclei 



L INTRODUCTION 

Astrophysical black holes span a large range of masses, 
from the remnants of stellar evolution to monsters weigh- 
ing by themselves almost as much as a dwarf galaxy. 
Notwithstanding the several orders of magnitude differ- 
ence between the smallest and the largest black hole 
known, all of them can be described by only two pa- 
rameters: mass and spin. So, besides their masses, M, 
astrophysical black holes are completely characterized 
by their dimensionless spin parameter, a = J hi J max = 
cJh/GM^, where is the angular momentum of the 
black hole, and < a < 1. 

Many theoretical efforts have focused on the mass 
growth of massive black holes (MB Hs) and on their 
feedback onto the host (see, e.g., IHaehnelt fc Ree^ 
1993; Haiman & Mcnou 2000; Kauffmann & Haehnelti 
~000£ .Wvithe fc Loeb. .2003^ Volonteri et al.i i20da 
Hopkins et all 120061 : iCroton et al.l 120061 ). Spin has 
received less attention in the cosmologic al context 
jModcrski et al. 1998; Volonteri et al. 2005; Lagos et al.l 
[20091: iFanidak is et al.l 120111: iBaraussd 120121) . but it has 
great relevance for the overall growth of MBHs, as 
follows. 

First, geometrically thin and optically thick accre- 
tion discs radiate with an efficiency, erad, which is al- 
most equal to the mass-to-energy conversion efficiency, 
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Crad ~ e = 1 - Eisco, whcrc Eisco = a/1 - 2/(3risco) 
is the specific energy of the gas particle (in rest mass- 
energy units) in the innermost stable {— marginally sta- 
ble) circular orbit (ISCO), and nsco is the radius of this 
orbit in GA4/c^ units. This radius and therefore radia- 
tion efficiency depend solely on the BH spin a. Maxi- 
mal efficiency (e ~ 0.42) is achievable by discs rotating 
around maximally spinning BHs; it drops to ~ 0.06 for 
non-spinning BHs, and to ~ 0.05 for maximally counter- 
rotating BHs. This entails a dependence of the BH mass- 
growth rate on the spin value, implying longest growth 
time scales for larger positive spins. More precisely, for 
a hole accreting at the Eddington rate, the black hole 
mass increases with time as: 

M(t) = M(0) exp {^—1-1-^ , (1) 

where igdd = M^^c^ / L^m = 5;^^ = 0-45 Gyr- The 
higher the spin, the higher e, implying longer timescales 
to grow the MBH mass by the same number of e-foldings. 

The radiative efficiency is also th e fundamenta l free 
parameter in the Soltan argument (iSoltanI I1982D and. 
more recently, in synthesis models (e.g.. lMerloni fc Heind 
l2008f ) which relate the integrated MBH mass density to 
the integrated emissivity of the AGN population, via the 
integral of the luminosity function of quasars. If the av- 
erage efficiency of converting accreted mass into luminos- 
ity is e = L/Mc^, then the MBH will increase its mass 
by M — {\ — e)Min, accounting for the fraction of the 
incoming mass that is radiated away. Applying this ar- 
gument to the whole MBH population, the MBH mass 
density, pbh, can be related to the integral of the lu- 
minosity function of quasar, ^'(L,2:), with the radiative 
efficiency being a free parameter: 
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Pbh{z)= —dz ^ ^^^{L,z)dL. (2) 

Jz "-^ Jo 

Recent results suggest that this approach might be 
too simplistic, as the radiative efficiency evolves along 
the cosmic time. iWang et al.l {2009) for instance sug- 
gests that quasars at the peak of their activity (z ^ 2) 
have high radiative efficiencies, hence high spins. At later 
times {z < 1), however, the average radiative efficiency 
decreases, hinting to lower spins. This paper addresses 
these issues, by providing spin distributions of MBHs as 
a function of cosmic epoch. 

MBH spins also affect the incidence of MBHs in galax- 
ies, via the "gravitational recoil" mechanism. When the 
members of a black hole binary coalesce, the center of 
mass of the coalescing system recoils due to the non- 
zero net linear momentum carried away by gravitational 
waves in the coalescence. If this recoil were sufficiently 
violent, the merged hole would break loose from the host 
and leave an empty nest. Recent breakthroughs in nu- 
merical relativity have allowed reliable computations of 
black hole mergers and recoil velocities, taking the effects 
of spin into account. Non-spinning MBHs, or binaries 
where MBH spins are aligned with the orbital angular 
momentum are expected to recoil with velocities below 
200 kms~^. The recoil is much larger, up to thousands 
kms~^, for MBHs with large spin s in non-aligned con- 
figurations (iCam panelli et al. 2067 HGonzalez et al.ll2007l : 
iHerrmann et al.ll2007( ). 

Finally, the spin of a hole might determine how 
much energy is extract able from the hole itself 
(jBlandford fc Znaiek|[l977l ) (Tchekhovskoy & McKinney 
2012). The so-called "spin paradigm" asserts that pow- 
erful relativistic jets are produ ced in AGN s with fast 
rotating black holes (Blandfor d et al.l I1990I ). implying 
that MBHs rotate slowly in radio-quiet quasars, which 
repre sent the majority of quasars ^Wilson fc Colbert! 
I1995D . Sikora et al. (2007) proposed a "spin-accretion 
paradigm" , suggesting that the production of powerful 
relativistic jets is conditioned by the presence of fast ro- 
tating holes, while it also depends on the accretion rate 
and on the presence of disc magneto-hydrodynamical 
winds required to provide the initial coUimation of the 
central Poynting flux dominated outflow, as in, e.g., the 
Blandford-Znajek process. 

As described above, MBH spins determine directly the 
mass-to-energy conversion efficiency of quasars. On the 
other hand, accretion determines the evolution of MBH 
spins. A hole that is initially non-rotating gets spun up 
to a maximally- rotating state {a — 1) after increasing 

its mass by a factor -v/G ~ 2.4. A maximally-rotating 
hole is spun down by retrograde accretion to a = after 
growing by a factor \/3/2 ~ 1.22. Different modes of 
MBH feeding imply different spin histories. Spin-up is 
a natural consequence of prolonged disc-mode accretion: 
any hole that increases substantially its mass by cap- 
turing material with constant angular momentum axis 
would ends up spinning rapidly ("coherent accretion"). 
Spin-down occurs when counter-rotating material is ac- 
creted, if the angular momentum of the accretion disc 
is strongly misaligned with respect to the direction of 
the MBH spin. It has been suggested that accretion 



may proceed also via small (i.e., the accreted mass is 
a very small fraction of the MBH mass, ^ 1% or less) 
and short uncorre l ated episodes ("chaotic acc retion" , 
IModerski fc Sikoral [T996t IKing fc Prin"gi^ [20061 ) . where 
accretion of co-rotating (causing spin-up) and counter- 
rotating (causing spin-down) is equally probable. As 
the IS CO for a retrograde orbit is at larger radii than 
for a prograde orbit, the transfer of angular momen- 
tum is more efficient in the former case. Accretion of 
counter-rotating material therefore spins MBHs down 
more efficiently th an co-rotating material spins them up. 
IKing et al.l (|2008D considered a MBH evolution scenario 
where chaotic accretion very rapidly adjusts the hole's 
spin parameter to average values a ~ 0.1 — 0.3 from a 
broad range of initial conditions, only weakly dependent 
on the overall angular momentum distribution of the ac- 
creting gas parcels. 

MBH-MBH mergers a lso influence the spin evolution. 
iBerti fc Volonteril ()2008l) consider how the dynamics of 
BH mergers influences the flnal spin. Except in the case 
of aligned mergers, a sequence of BH mergers can lead 
to large spins > 0.9 only if MBHs start already with 
large spins and they do not experience many major merg- 
ers. Therefore, the common assumption that mergers be- 
tween MBHs of similar mass always lead to large spins 
needs to be revised. 

The focus of this paper will be on the cosmic evolu- 
tion of spins of massiv e black holes, M ^ 10^ — 10^ Mq 
(|Richstone et al.l 119981: iFerrarese fc Ford"2005V specifi- 
cally on how accretion and MBH-MBH mergers deter- 
mine the magnitude of spins. In section 2 we describe 
the basic infrastructure that we use to model the cosmic 
evolution of structures. In section 3 we summarize how 
we model spin evolution in MBH-MBH mergers, while 
in section 4 we describe how different phases of accre- 
tion, related to the cosmic evolution of galaxies and of 
the MBHs they host, influence MBH spins. In section 5 
we consider a series of observational constraints that we 
adopt to anchor our model. In section 6 we discuss our 
results, and we present our conclusions in section 7. 

2. THE BACKBONE: DARK MATTER HALOS AND 
GALAXIES 

We investigate the evolution of MBHs via cosmological 
realizations of the merger hierarchy of dark matter halos 
from early times to the present in a ACDM cosmology 
(WMAP5, Komatsu et al. 2009). We track the merger 
history of 300 parent halos with present-day masses in 
the range 10" < Mh < 10^'' with a Monte Carlo 
algorithm (I Volonteri et al.ll2003D . The mass resolution 
of our algorithm reaches IO^Mq at z = 20, and the most 
massive halos are split into up to 600,000 progenitors. 

We wish to keep our models as simple as possible, 
while making sure that the properties of the MBHs we 
study are correctly determined through their cosmic evo- 
lution of their hosts. We do not explicitly model the 
evolution of the baryonic component of the host galaxies 
through cooling, star formation and va rious feedbacks 
fsee jLagos et al.l 120091 : [Fanid akis et al.l 12011: Baraussi 
120121 and references therein for models that treat in de- 
tail semi-analytically the baryonic component of galaxies 
and its link to MBH evolution). In our models we use 
only one parameter to link the host halo to the central 
MBH, and it is the host's central velocity dispersion. We 
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link the central stellar velocity dispersion of the host to 
the asymptotic virial velocity (wvir) assuming a spherical, 
isothermal halo, so that a = Vvir//^- We calculate the 
circular velocity from the mass of the host halo and its 
redshift. A halo of mass Mh collapsing at redshift z has 
a circular velocity: 
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where Ac is the over-density at virialization relative 
to the critical density. For a WMAP5 cosmology we 
adopt here the fitt ing formula Ac = IStt^ + 82d — SOrf^ 
(|Brvan fc Norman|[l998 ). where d = - 1, and f2,^ = 
^m(X + zfl{^^(\ + zf + + r!fe(l + zf). 

At high redshift we seed dark matter halos with MBHs 
created by gas collapse. Specifically, we adopt here 
the formation model detailed in Natarajan & Volonteri 
(2011) based on Toomre instabilities (Lodato & Natara- 
jan 2006) . The Toomre parameter is defined as Q = , 
where S is the surface mass density, Cs is the sound speed, 
K — \/2V\i/R is the epicyclic frequency, and V is the cir- 
cular velocity of the disc. When Q approaches a critical 
value, Qc, of order unity, the disc is subject to gravita- 
tional instabilities. If the destabilization of the system is 
not too violent, instabilities lead to mass infall instead of 
fragmentation into bound clumps and global star forma- 
tion in the entire disk (Lodato & Natarajan 2006). This 
process stops when the amount of mass transported to 
the center is enough to make the disc marginally sta- 
ble. The mass that has to be accumulated in the center 
to make the disc stable, Minf, is obtained by requiring 
that Q = Qc- This condition can be computed from 
the Toomre stability criterion and from the disc prop- 
erties, determined from the dark matter halo mass, via 

Tvir oc Af^^^, and angular momentum, via the spin pa- 
rameter, AsDin: 
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for A < A„,ax = /dQc/8(/d/id)(7;ir/Tgas)i/^^ Here A,„ax 
is the maximum halo spin parameter for which the disc 
is gravitationally unstable, fd = 0.05 is the gas frac- 
tion that participates in the infall, and jd — 0.05 is the 
fraction of the halo angular momentum retained by the 



collapsing gas. We further assume Tg, 



5000K and 



Qc = 2 (see Volonteri, Lodato & Natarajan 2008 for a 
validation of the parameter choice). Given the mass and 
spin parameter of a halo, the mass that accretes to the 
center in order to make the disc stable, Minf., is an upper 
limit to the mass that can go into MBH formation. We 
assume here Mgeed = Minf. Please refer to Natarajan & 
Volonteri (2011) and references therein for details on the 
MBH formation process. 

To model further MBH growth, for simplicity we fur- 
ther assume that the black hole-tr (M-cr) scaling is: 



= (200^^ 
(jTremaine et al.|[200l . 
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Our model for MBH growth requires only to know 
whether a galaxy is gas-rich (we typically refer to 
gas-rich galaxies as "discs" in this paper) or gas-poor 
("spheroid"). Since forming galaxy discs even in high- 
resolution cosmological simulations is extremely chal- 
lenging (iGovernato et al., 2007i ). we keep our model for 
galaxy morphology as simple as possible. Morphology 
is related to the merger history, using a three-parameter 
model, where spheroid formation depends on both halo 
mass ratio and the absolute halo mass, and a spheroid 
can re-acquire a disc through cold flows and mergers 
with gas-rich galaxies. iKoda et al] ()2007D show that 
the fraction of disc- vs spheroid-dominated galaxies is 
well explained if the only merger events that lead to 
spheroid formation have mass ratio >0.3 and virial ve- 
locity > 55kms~i; also, the merger timescale must be 
less than the time between when the merger starts and 
today, z = 0. We assume that spheroids form after a 
merger that meets these requirement. We additionally 
allow a disc to reform after 5 Gyrs in galaxies with virial 
velocity < 300 km s^^ where no major mergers occurred 
to include the effect of cold flows. In section 5 we discuss 
the sense of this approach. 

3. SPIN EVOLUTION DUE TO MBH-MBH 
MERGERS 

We assume that, when two galaxies hosting MBHs 
merge, the MBHs themselves merge within the 
merger tirnescale of the host halos (jSesana et al.l 120071 : 
iDotti et al .''2007'. and references ther ein). We adopt the 
relations suggested by Bovlan-Kolch in et al.l (j2008f ) for 
the merger timescale. We model MBH spin changes due 
to mergers adopting an analytical scheme similar to that 
described in Berti & Volonteri (2008), based on simu- 
l ations of black hole mergers in full ge n eral relativity 
(iRezzolla et al.l 120081: ITousto et al.ll2009[ ). iKesden et atl 
( 2010() has validated the consistency of different fitting 
formulae for calculatin g the spin of MBH r emnants, and 
we refer the reader to ILousto et al.l ()2010[ ) for the most 
comprehensive fitting formulae. Due to computational 
constraints, we adopt here the fitting formulae of Rez- 
zolla et al. (2008) for their easy implementation: 
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2 (|ai| cos/3 + |a2|q2 C0S7) \i\q + ,(6) 

where q = M2 /Mi ^ 1 is the mass ratio between the two 
MBHs; the three (cosine) angles a, (3 and 7 are defined 
by 

cos a = Oi • 02 , cos (3 = ai ■ I . cos 7 = 0,2 • € . (7) 

where |£| is the magnitude of the orbital angular momen- 
tum, and 



1^1 = 



S4 



(l + g2)2 
s^v -I- to + 2 



(|ai|2 + |a2p/ + 2|ai||a2|g^cosa) + 



l + g2 



(|ai| cos/3 -I- |a2|g^ cos 7) 



(8) 
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where v is the symmetric mass ratio v = M1M2/ {Mi + 
M2)^, and the coefficients take the values S4 — — 0.129 ± 
0.012, S5 = -0.384 ± 0.261, to = -2.686 ± 0.065, ^2 = 
-3.454 ± 0.132, h = 2.353 ± 0.548. 

We model the spin-orbit configuration differently de- 
pending on the properties of the host galaxies. Dur- 
ing a gas-rich merger, large amounts of gas are driven 
toward the centers of t he two interacting galaxies 
(jPownes fc Solomonlll998h and form a dense circumnu- 
clear disc in which the MBHs settle. In this phase, 
the MBHs accrete in a coherent manner at a rate suf- 
ficient to align their spins, initially oriented at ran- 
dom, to the angular mo mentum of the nuclear disc (jLiul 
120041: IBogdanovic et al.ll2007t IDotti et al.ll2010D : in re- 
sponse to the Bardeen-Petterson warping of the small- 
scale accretion discs grown around each MBH, total an- 
gular momentum conservation imposes fast (<1 Myr) 
alignment of the BH spins with the angular momen- 
tum of their orbit and so of the large-scale circumnu- 
clear disc (Dotti et al. 2010). Thereafter, accretion re- 
mains prograde until coalescence, with no major changes 
in the MBH spin orientation. Under these circumstances, 
cos a = cos /3 = cos 7 = 1, and the MBH remnant retains 
the spin direction of the parent MBH spins, both ori- 
ented parallel to the angular momentum of their orbit: 
the post-coalescence MBH may thus acquire a large spin 
> 0.7 - 0.9 (Berti & Volonteri 2008), sum of the internal 
and orbital spins. 

In the case of gas-poor mergers, instead, i.e. when 
MBHs do not accrete during mergers, evolving solely via 
stellar dynamical processes, there is no reason to exp ect 
any symmetry or alignment (jBogdanovic et aDl2007| ). so 
isotropy should be a good assumption in the absence of 
accretion discs or gas, so that cos a, cos/?, a nd cos 7 are 
isotropically distributed. iBerti fc Volonteril (j2008il show 
that for isotropic configurations m ergers tend to "spin- 
down " a fast-spinning hole (see also lHughes fc BlandfordI 
120031 ). For intermediate-large mass ratios (mass ratio q = 
-^BH,2/-^^BH,2 ^ 1 between 0.1 and 1) mergers tend to 
produce MBHs with average spins very close to the value 
~ 0.7 resulting from equal-mass, non-spinning mergers. 
For smaller mass ratios the larger MBH dominates the 
dynamics, and the final spin can be substantially larger 
or smaller than this value. 

4. SPIN EVOLUTION DUE TO ACCRETION 

We discuss here the feeding of MBHs in the quasar 
phase and its aftermath and in the more quiescent Seyfert 
galaxi es. Simulations of galaxy mergers and MB H ac- 
tivity ()Di Matteo et al.ll2005t iHopkins et aIll2006D show 
that for every accretion episode triggered by a galaxy 
merger a three phase picture can be drawn. At the be- 
ginning the MBH has an "healthy diet", with /Edd = 
L/L^dA s% 1- When the MBH mass reaches the "M-cr" 
relation, the MBHs feedback can be sufficient to unbind 
and "blowout" the gas feeding it (jHopkins fc Hernquistl 
[2006I) . causing a final "starvation", when /Edd rapidly 
decreases, until no more gas is available to feed the 
MBH. We argue that during the healthy diet and blowout 
phases MBHs g ain a high spin wh ile accreting efficiently 
and coherently (jDotti et al.ll20ldl ). building the popula- 
tion of high-z quasars. Coherent accretion ensues because 
MBHs in merger remnants are expected to be surrounded 
by dense circum- nuclear discs (Sanders fc MirabcL1996 ). 



Maio et al. (2012) study the evolution of the angular 
momentum of material feeding MBHs embedded in cir- 
cumnuclear discs, and they find that coherence of the 
accretion fiow near each MBH reflects the large-scale co- 
herence of the disc's rotation. 

After the "blowout" phase, starving MBHs are no 
longer surrounded by a thick gas disc, that determines 
the angular momentum of the material ending up in the 
accretion disc. During this last phase we do not expect 
accretion to necessarily proceed coherently any longer. 
During the starvation phase the accretion rate decreases 
rapidly. The depletion of gas in galaxies and the decrease 
in the galaxy interaction rate at late cosmic times causes 
therefore a widespread "famine" in low- 2; ellipticals (i.e., 
the merger remnants), where AGN with low accretion 
rates dominate ("radio" m ode, see, e.g., ICroton et al.l 
f200l lChurazov et al.l [2005'). 

We investigate the evolution of MBH s pins during the 
quasar pha se expanding previous work (jVolonteri et al.l 
120051 [2007|) to more realistic models. We model the joint 
mass and spin evolution by coupling the results on the 
mass accretion rate as a function of time in simulations 
(Hopkins et al. 2005; Volonteri et al. 2006,) , with the spi n 
evolution due to disc accretion (|Volonteri et al.ll2007D . 
thus solving a system of two coupled differential equa- 
tions (/Edd as a function of M and time; spin a as a 
function of /Edd, M . The framework has been derived in 
Volonteri et al. 2005, Volonteri et al. 2006, Hopkins fc 
Hernquist 2006, and Volonteri et al. 2007). We remind 
here the relevant information. 

4.1. Quasar phase 

After a halo merger with mass ratio larger that 3:10, 
in which at least one of the two is a disc galaxy (hence, 
with conspicuous cold gas content) we assume that a 
merger-driven accretion episode is triggerecj^. After a 
dynamical timescale is elapsed (roughly, after the first 
pericentric passage, cf. I Van Wassenhove et al.l [20121 ) ac- 
cretion starts. If at that point the MBH mass lies below 
the M-cr relation, accretion occurs at the Eddington rate 
(/Edd = 1)- Additionally, during this early phase the 
MBH is nested into a nuclear disc that feeds the MBH 
coherently. The following scheme is applied to the joint 
evolution of mass, spin and radiative efficiency in this 
phase. Let us define M and a as the black hole mass and 
spin parameter at the beginning of the timestep, and /i 
as the cosine of the angle between the MBH spin and 
the inner accretion disc angular momentum. Irrespec- 
tive of the infalling material's original angular momen- 
tum vector, Lense-Thirring precession imposes axisym- 
metry close in, with the gas accreting on either prograde 
(/i = 1) or retrograde equatorial orbits (/i = — 1). In 
natural units, where c — G — 1: 



nsco = 3 + Z2 T v/(3-^i)(3 + Zi+2Z2), (9) 

^ Please note that this assumption is at variance with pre- 
vious models of MBH cosmic evolution within ou r framework 
(e.g., Volonteri ct al. 2003, 2008; Volonteri & Nata raianI |2009| : 
Volonteri ct al. 2009; Volonteri & Be gelman _2010), where the 
threshold was set to a lower value of 1:10. The reason for in- 
creasing the threshold for QSO phase to occur is the addition of 
avenues for MBH growth other than merger-driven accretion and 
MBH-MBH mergers. 
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is the radius of the ISCO, where Z\ and Z2 are functions 
of a only (Bardeen, Press, & Teukolsky 1972), 

= 1 + (1 - a2)i/3[(l + afl^ + (1 - af'X (10) 
Z2 = [3a2 + z2]i/2^ (11) 

and the upper (lower) sign refers to prograde (retrograde) 
orbits. We calculate the accretion efficiency as: 



e = 1 — -E'isco, 
2 



ISCO 



= 1 



Srisco 



1/2 



(12) 
(13) 



which is is also a plausible assumption for the radia- 
tive efficiency (e, mass-to-energy conversion) for thin- 
disc accretion occurring at large fractions of the Edding- 
ton rate (see below for the case of radiatively inefficient 
flows). We calculate self-consistently the radiative effi- 
ciency from the MBH spin and the location of the ISCO 
(i.e., taking into consideration the direction of the rela- 
tive angular momentum of spin and disc, co- or counter- 
rotating) . 

Assuming that during time a timestep Ai ^ 10^ — 10^ 
yr the radiative efficiency and Eddington rate remain 
constant, e=const and /Edd=const, from the derivation 
shown in the Appendix, one obtains that the MBH mass 



grows as: 



M{t + td) = M{t) cxp /j 
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where tEdd 



0.45 Gyr and /Edd represents 
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the Eddington fraction. We update the magnitude of 
the MBH spin through: 
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M{t) 
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a{t + At) - 0.998 for > rll^oi*) (16) 

(Bardeen 1970 After updating the spin magnitude, we 
also update the mass-to-energy conversion efficiency by 
determining the new ISCO corresponding to a{t + At), 
and therefore e{t + At) to be used at the successive 
timestep iteratively. We here assur ne fast alignment be- 
twccn accretio n disc and s pin (see 'Nataraian & Pringle 
[1998 : Volontc ri et al.|[2005l . 12007; .Pere go et a l. 2009), as 

We limit the MBH spin to a = 0.998 following the calculation 
of Thorne (1974) that showed that the radiation emitted by the disc 
and swallowed by the hole produces a counteracting torque, which 
prevents spin up beyond . We note that magnetic fields connecting 
material in the disc and the plunging region may further reduce 
the equilibrium spin by transporting angular momentum outward 
in non-geometrically thin discs. Fully relativistic magnetohydrody- 
namic simulations for a series of thick accretion disc models show 
that spin equilibrium is reached at a Ri 0.93 (Gammie et al. 2004), 
while in slim discs acc retion can increase the spin up to a = 0.9994 
aSadowski et al.|[20TlD 



the alignment timescale is ~ Myr, so that accretion is 
prograde during the quasar phase. 

4.2. Decline phase 

When a MBH reaches a mass close to the value 
corresponding to the M-a co rrelation (Mbh .<j) for 
its host, we assume , follo wing IHopkins et al.l (|2006); 
iHopkins fc Hernguistl ()2006D that a self-regulation ensues 
and the MBH feedback unbinds the gas closest to the 
MBH, thus reducing its feeding. To match the luminos- 
ity function of quasars, we start the decline phase when 
the MBH mass is 0.25 x MsH^a- The MBH continues 
accreting during the decline of accretion and it typically 
reaches MBH,a by the end of the accretion episode. 

We model the decreas e of the accretion rate throu gh 
the analytical formula in IHopkins fc Hernguistl (|2006[) : 



fEdd{t) 



t + t 



t 



/e 



(17) 



where and tf^^^ ~ 4.1 x lQ^{MBH^a /lG^)y^e yr, 

and t — Q (where /Edd — 1) represents the time when 
the MBH reaches the threshold (0.25 x MBH,a)- Since 

Mia{t) — fEddit) M{t)/tEdd , depends on time, the ac- 
cretion rate must be integrated self-consistently, and the 
mass now grows with time as: 



M(t) = M(0)exp 




1 



(18) 



We use again Eq. [91IT61 to model spin evolution, how- 
ever, we explore two possible scenarios. In our reference 
case we assume that the "outflow" causes some stirring 
of the angular momentum of the gas within the central 
region. We therefore explore a "chaotic" case where in 
the decline phase we pick a new random /i = lor/i = — 1 
at each timestep (~ 10^ — 10^ years) to mimic the lack of 
coherence in the accretion flow after MBH feedback has 
blown away the surrounding gas. In a second model we 
assume instead that the "blow-out" of gas does not af- 
fect strongly the angular momentum of the material near 
the MBH, and persist with keeping /i = 1. Given that 
the timescale for decline is longer for larger MBHs (cf. 
Eq. ITTI) , the larger the MBH the longer the phase at rel- 
atively high accretion rates, and the faster the decrease 
of MBH spin, as more mass is accreted in a non-coherent 
fashion. 

When the accretion rates become very sub-Eddington, 
we assume that the accretion flow becomes optically 
thin and geometrically thick. In th is state the radiative 
power is strongly supp ressed (e.g.. iNaravan fc Yil 119941 : 
lAbramowicz et al] 1 19951 ). so that the radiative efficiency 
differs from the mass-to-energy conversion efficiency, e, 
that depends on the location of the ISCO only. Indeed, 
the radiative efficiency becomes very model dependent 
and uncertain. In order to estimate the effect of radia- 
tively inefficient accretion on the MBH population we 
adopt here a speci fic functional form for th e radiative ef- 
ficiency. Following iMerloni fc Hein3 ()2008l ) we write the 
radiative efficiency, eiad, as a combination of the mass-to- 
energy conversion, e, and of a term that depends on the 



6 



Volonteri et al. 



prope rties of the accretion flow itself. iMerloni fc Heinj 
(|2008D suggest that the transition in the disc proper- 
ties occurf@ at /Edd < /sddxr = 3 X 10~^, and that 
Crad = e (/Edd//Edd,cr)- This spccific choice allows us to 
estimate qualitatively the impact of radiatively inefficient 
sources to the AGN populations, and on the inferences 
that one can (or not) make from observables. 

4.3. Quiescent elliptical phase 

After the formation of an elliptical galaxy, the feed- 
ing of the MBH can be sustained by the recycled gas 
(primarily from red giant winds an d planetary nebu- 
lae) of the evolving stel lar population jCiotti fc Ostrikeil 
■1997, 2001, 2007; Ci otti et al.l [20T0I ). As shown by 
Ciotti fc Ostrikeii (120 llT the behaviour of the accretion 



rate is similar to what we describe above: at early times 
the evolution is characterized by major, albeit intermit- 
tent, accretion episodes, while at low redshift accretion is 
smooth and characterized by /Edd ^ 1- We have imple- 
mented this channel of MBH feeding only for quiescent 
ellipticals, that is only after a spheroid is formed and had 
time to relax. When the MBH mass is ^ 10~^ of the 
spheroid mass (iMagorrian et al.l 119981: iMarconi fc Huntl 
120031 iHaring fc Rixl I2004D . a'^ the spheroid is mod- 
eled as a Hernquist p rofile, the geometrical model by 
IVolonteri et al.l (j2011a[ ) implies that the quiescent level 
of accretion onto a central MBH due to recycled gas is 
/Edd — 10~^. We assume here that this mode of accretion 
has a constant /Edd = 10~^, and we model spin evolution 
assuming that that the accreted material is isotropically 
distributed (i.e., we pick a random /i = 1 or /i = — 1 at 
each timestep, and use Eq. to evolve the spin mag- 
nitude). This phase has little effect on the magnitude 
of MBH spins, due to the low accretion rates and the 
modest implied MBH growth (cf. Eq. [T| 



4.4. Molecular cloud accretion in disc galaxies 

Several observations suggest that single accretion 
events last ~ 10^ years in Seyfert galaxies, while the to- 
tal activity lifetime (based on the fraction of disc galax- 
ies that are Seyfert) is 10* — 10^ years (e.g.. lKharb et al.l 
120061 : IHo et al.llT997D . This suggests that accretion events 
are very small and very 'compact'. Smaller MBHs, 
powering low luminosity AGN, likely grow by accret- 
ing smaller packets of mater ial, such as molecular clouds 
(|Hopkins &: Hernquisti 120061 ). Compact self-gravitating 
cores of molecular clouds (MC) can occasionally reach 
subparsec regions. In gas-rich, star-forming disc galaxies 
the MBH is likely to be fed by short, recurrent, uncorre- 
lated accretion episodes. The spin evolution of a MBH 
hosted by a quiescent disc galaxy would then resemble 
the "chaotic accretion" scenari o. This argum e nt wa s 
discussed only qualitatively by IVolonteri et al.l ()2007[ ). 
We now wish to provide quantitative statistical predic- 
tions for the di stribution of M BH s pins in different hosts. 
We fo llow here lSandersI (jl98lD and lHopkins fc Hernquisti 
(|2006D to determine the event rate, and Volonteri et al. 
2007 to couple accretion episodes to spin evolution. 



8 Please note that IMerloni k. Heind ll2008t) use a different no- 
tation and terminology. Their A is our /Edd = ^/^Edd = 
Edd: and their fn = e/Edd/^rad- As long as the ac- 
cretion flow is optically thick and geometrically thin, i.e., before 
the transition to very sub-Eddington flows, m = /Edd- 



In a disc galaxy, at each timestep. At, we determine 
the probability of a MG accretion event as: 



At lO^^g.^^ 



tMC 



Rc 



(19) 



where Rc\ 10 p c (Hopkins & Hernquist 2006). As 
in Vo lonteri et al.l (2007) we further assume a lognor- 
mal distribution (peaked at log(AfMc/Mg) = 4, with 
a dispersion of 0.75) for the mass function of MG close 
to galaxy centers (based on the Milky Way case, e.g., 
IPerets et al] [20071 ). 

We model accretion of MGs through a description in- 
spired by Sanders & Bottema (1983) and Wardle & 
Yusef-Zadeh (2008). We assume that the MBH cap- 
tures only material passing within the Bondi radius, Rb, 
and we also assume that specific angular momentum is 
conserved, so that the outer edge of the disc that forms 
around the MBH corresponds to the material originally 
at the Bondi radius: 



Rd = 2X^Rb = 8.9 pcA^ 



M 



lO^Mo VlOOkms-i 



, (20) 



where A is the fraction of angular momentum retained 
by gas during circularization. The maximum captured 
mass will be contained in a cylinder with radius Rb and 
length 2 x i?ci, the MG diameter. If k is the ratio of the 
mass going into the disc with respect to the whole mass 
in the cylinder, then: 



^k^Mmc =4.7x 10^ 
Rcl 



f ^ 

VlOOkms-1 



M^ 



(21) 



The inflow time will be of order of the viscous timescale 
for the disc. 



alGMd 
a 



1/2 



(lOOkms-i) 



3.4 X 10^ 



,A^_M_ 



(22) 



-lOOkms-i. 

so that for the whole disc to be consumed we can cal- 
culate an upper limit to the mean accretion rate and 
luminosity: 



Mmax = —r =0.13- 



A3 lO^M^ 



(1001^)"^-^-^^"'' 



(23) 



where iV23 is the column density in the MG in units of 
10^'^ cm~^, and a.v — 0.1 is the viscosity parameter. Fol- 
lowing Wardle & Yusef-Zadeh (2011) we set A = 0.3 and 
K = 1. If Mmax is less than the Eddington rate (assuming 
a radiative efficiency of 10%) we let the MBH accrete the 
whole M(j max over a time tviscj otherwise we treat accre- 
tion similarly to the "decline" phase of quasars (Eq. ITTI 
andll8p. as feedback from the high luminosity produced 
by accreting the cloud will limit the amount of material 
the MBH can effectively swallow. 

From Eq.[2Tl it is evident that the mass accreted in one 
of these episodes is typically much less than the mass of 
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Figure 1. Properties of all MBHs fed by MCs at all redshifts: 
normalized distribution of the masses of MBHs (top) and their 
host halos (bottom). The probability of MC accretion increases 
with galaxy mass, on the other hand most massive galaxies are 
spheroids, and therefore they do not have a population of MCs 
available. 

the MBH (typically between lO'^ and 10"^ of the MBH 
mass), we therefore assume that no alignment between 
accretion disc and MBH spin can occur, and that retro- 
and prograde accretion is equally probable, i.e., we as- 
sign n — 1 or /i = —1 with equal probability at each 
event, and keep fi constant over the accretion phase, us- 
ing Eq. [16] to evolve the spin magnitude. As shown by 
Volonteri et al. 2007 the os these assumptions result is a 
spin down in a random walk fashion that depends on the 
mass of the MBH and on the number of events. We note 
that this is an extreme condition of randomness in the 
orbits and distribution of MCs, as any common sense of 
rotation caused by the presence of disc-like structures in 
the host would decrease the degree of anisotropy (Dotti 
et al. 2012). 

Eq. [m instead, shows that accretion of MCs is more 
probable in large galaxies, since the accretion probabil- 
ity is directly proportional to the velocity dispersion of 
the galaxy. Large galaxies, however, are more likely to 
be gas-poor spheroids. Large galaxies therefore have a 
lower probability of being gas-rich and host a population 
of MCs, while at the same time their MBHs have a higher 
probability of capturing MCs, if clouds are present. 
Thus, this type of accretion events occur typically in 
galaxies hosted in halos with mass ~ 10"'^^ — 10^^ Mg and 
fuel MBHs with mass 10^ - 10^ Mq. In Fig. [1] we show 
the distribution of the masses of MBHs and their host 
halos where accretion of MCs takes place according to 
our scheme. 

5. ANCHORING THE MODEL 

Here we present the constraints that we use to "an- 
chor" our model to observed properties of galaxies and 
AGN. After validating our scheme for accretion and host 



evolution, we will discuss what its implications are for 
the, still unknown, distribution of MBH spins. 

First, in Fig. [2] we compare our model to constraints 
at z = 0. In the bottom panel we show the morpho- 
logical fraction as a function of galaxy stellar mass. We 
scale from halo m a.ss to stellar mass thr ough the data de- 
scribed in Fig. 1 of lHopkins et al.l ()2010[ ). assuming that a 
fixed fraction 10-50% of the baryons is in stars (open red 
points and filled orange ones respectively). While agree- 
ment is far from perfect, our simple approach reproduces 
the correct trend. We recall here that we do not model 
the whole evolution of gas and stars, nor disc formation 
(cf. lBaraussil2012l for a comprehensive model), but to 
derive all properties our scheme uses only a single quan- 
tity, the halo mass. 

In the top panel of Fig. [2] we reproduce the re- 
lationship between MB H masses and circular velocity 
(jVolontcr i et al.l |2011bl). Circles are model MBHs at 
z = 0, while errorbars show the data (Kormendy & 
Bender 2011; Volonteri et al. 2011). We note that 
our model fails to reproduce some of the massive MBHs 
at = in spheroids. The reason for the suppressed 
growth is that no MC accretion can happen in our scheme 
in gas-poor galaxies, and we have not implemented a 
time-dependent a. c cretio n of recycled gas. As noted by 
iCiotti fc Ostrike^ pOOTi ) the behavior of a MBH fueled 
by stellar mass loss is self-regulated between "on" and 
"off' phases. Our model includes only the "ofP' (quies- 
cent) phase and therefore underestimates the growth due 
to this fueling channel. 

Our second anchor is the luminosity function of AGN 
in the redshift range 0.5 ^ z ^ 3. This is a strong con- 
straint for our accretion scheme, and it is shown in Fig.[3l 
Our scheme produces an AGN population in good agree- 
ment with the observations at z — 0.5, z — 2 and z = 3, 
while we slightly underproduce AGN at z = 1. The lack 
of high luminosity quasars at z = 2 and z = 3 is due 
to our merger trees not including halos massive enough 
to host MBHs with masses above a few 10* Mgat those 
redshifts. Overall, however, we obtain the correct trends. 
This figure shows that there is little difference in the two 
models we explore on the effect of feedback over the an- 
gular momentum of the nuclear gas ("chaotic" or "co- 
herent" decline phases). The difference between chaotic 
and coherent decline (the "quasar" phase is coherent in 
all cases), reflects only on spin and as a consequence on 
radiative efficiency, not on the Eddington ratio. In the 
following we will distinguish the two models only when 
the results are significantly different, otherwise we show 
only the reference case ( "chaotic" decline phase) . 

The main parameters influencing the performance of 
the model against the constraints are the mass ratio 
above which a mergers can trigger quasar activity and 
the fraction of MBH,a when the decline phase starts. The 
former is set to > 3 : 10 in the present model to match as 
well as possible the luminosity function at the bright end 
without overestimating MBH masses at a given circular 
velocity at z = 0. We note that MC accretion accounts 
for almost all sources up to L ~ 10^^ L© aX z — 0.5 — 1 
and L ~ 10^° L© at z = 2 — 3. For the latter parameter, 
we adopted a value of 0.25 x MBH,a to avoid overesti- 
mating the MBH masses at z = 0. 

Finally, we generate ad hoc merger trees of Mh = 
2 X 10^^ Mq halos at z = 5 — 6 to check that our 
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Figure 2. Top panel: relationship between MBH masses and 
circular velocity. Circles are model MBHs at z = 0, errorbars 
show the datapoints co llected in Kormendy & Bender (2011) and 
I Volonteri et al .l ((201TB!). Bottom panel: fraction of spheroids as a 
function of stellar mass. We scale from halo mass to stellar mass 
through the data described in Fig. 1 of Hopkinset al. 2010, assum- 
ing that a fixed fraction 10-50% of the baryons is in stars (open 
red points and filled orange ones respectively). We compare the 
fraction of spheroids to Conselice et al. (2006). 

model reproduces the existence of powerful quasars at 
z ^ 6, and that the mass density we obtain at z > 5, 
of 6 X 10"^ Mg/Mpc"^, does not overproduce the X-ray 
background (upper limit of 10^ Mg/Mpc'^, Salvaterra et 
al. 2012). Fig. |4] compares the theoretical mass function 
of MBHs that power quasars with bolometric luminosity 
larger than 10^^ er g/s at z = 6 t o the empirical mass 
function derived bv lWillott et al.l ()201Clf ) from a sample 
of z = 6 quasars in the Canada-France High-z Quasar 
Survey. 

As discussed above, our model, while far from being 
able to explain every single detail of the MBH population 
and its growth, qualitatively grasps most of the global 
behavior. We therefore consider our attempts to model 
the spin evolution also of qualitative nature. Regardless 
of the simplified nature of our models, we can learn how 
different patterns influence the evolution of MBH spins. 

6. ACCRETION AND SPIN EVOLUTION: RESULTS 

In Figure [5] we show examples of spin evolution of the 
MBH hosted in a large spheroid today along its cosmic 
history. Most of the accreted MBH mass is accumulated 
during episodes of efficient growth at early times (up un- 
til z ~ 2). At lower redshift the MBH grows mostly 
through MBH-MBH mergers. While different prescrip- 
tion for "chaotic" or "coherent" post-feedback phases 
lead to different histories for the MBH spin, the final 
spin is set by a MBH merger between two roughly equal 
mass systems occurring at z ~ 0.3. In the case of disc 
galaxies, late phases of MC accretion at substantial ac- 
cretion rates (> 10~^ in Eddington units ) contribute to 
setting the final spin of the MBH (Fig. 
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Figure 3. Luminosity function at different redshifts. We show 
here minimum and maximum values, considering both l-<j sta- 
tistical uncertainties, using Poissonian statistics, and fraction of 
absorbed AGN (La Franca et al. 2005). Orange (45° hatching): 
coherent accretion during the decline of the quasar phase. Gray 
(horizontal hatching): chaotic accretion during the decline of the 
quasar phase. 
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Figure 4. Red histogram: theoretical mass function of MBHs at 
z = 6 that power quasars with bolometric luminosity larger than 
10**^ erg/s. Blue points: mass function derived by Willott et al. 
2010 from a sample of z = 6 quasars. 



Statistically, Fig. [7] shows the evolution of the loga- 
rithmic Eddington ratio as a function of redshift. The 
top panel shows a sample selected on the MBH mass 
{M > 10'^Mq and M > lO^Mg) showing classic "anti- 
hierarchical" behavior (Merloni 2004) , with the the most 
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redshift 

Figure 5. Bottom panel; mass growth of a MBH in a galaxy that 
becomes a large spheroid by 2 = 0. The evolution is extracted 
from a merger tree describing a dark matter halo of mass 4 X 
10^"^ Mq at 2 = 0. The galaxy is the central galaxy in a group-sized 
halo. Middle panel: evolution of the Eddington rate vs redshift. 
Top panel; evolution of the spin parameter. After an early phase 
of rapid accretion and growth the accretion rate declines and the 
MBH is fed by stellar winds only (quiescent phase) in the past ~ 6 
Gyr. Today's spin is defined by a MBH-MBH merger, with a mass 
ratio q = 0.35 that occurred ~3 Gyr ago. Here the spin of the 
MBH before a MBH-MBH coalescence is shown as a blue empty 
circle, an the spin after coalescence as a red filled dot. 
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Figure 6. Evolution of a MBH in a galaxy that becomes a Milky 
Way- type disc by z = 0. The evolution is extracted from a merger 
tree describing a dark matter halo of mass 10^2 Mq at 2 = 0. The 
galaxy is the central galaxy in a Local Group-sized halo. Note the 
occurrence of MC accretion at z < 0.5. Panels, lines and symbols 
as in Fig. |5] 



massive MBHs being more active at earlier cosmic times. 
At a given luminosity threshold the typical Eddington 
ratio decreases slightly at late times (bottom panel of 
Fig. [T]) , tracking instead the overall increase in the mass 
of MBHs. 

In Fig. [8] we compare the spins of MBHs hosted in discs 
and spheroids, in different redshifts bins. At low MBH 
mass (M ~ 10^ Mg) MBHs hosted in gas-rich galaxies 
tend to have low spins. The spin distribution in gas-rich 
galaxies tends to move towards higher spins as mass in- 
creases. This is mostly related to MC accretion. For the 
most massive MBHs most of the growth occurs earlier 
on through merger-driven accretion that tends to spin- 
up MBHs. Accretion of MCs does not modify much the 
spins of these MBHs because the total angular momen- 
tum accreted through MCs is less than the total angular 
momentum the MBH has (i.e., the total mass accreted 
by the MBH through MCs is much less than the mass 
of the MBH). On the other hand, for low- mass MBHs 
the mass accreted in MCs is of the same order as the 
MBH mass, therefore they have a stronger effect on the 
spin distribution, lowering the typical spin of low-mass 
MBHs (we remind that we have assumed that MCs ac- 
crete isotropically on MBHs) . The distribution of spins of 
MBHs hosted in gas-poor galaxies has little dependence 
on mass and redshift. In these galaxies, in general, most 
MBHs have spin a ~ 0.4 — 0.8. Spins tend to slightly 
decrease as MBH mass increases. We find no strong de- 
pendence on whether accretion occurs mostly chaotically 
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Figure 7. Mean values of the Eddington ratio (logarithmic units) 
as a function of redshift. The top panel shows a sample selected 
on the MBH mass (lO'^Mo < Af < lO^Mg and M > IO^Mq) 
showing classic "anti-hierarchical" behavior. The bottom panel 
focuses instead on luminosity-selected AGN, for which the typical 
accretion rate is much higher. 
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or coherently after 'feedback' effects take place, except 
at the highest masses. We have run a test case where 
we have artificially "turned-ofF' spin evolution via MBH- 
MBH mergers (while keeping the mass increase through 
mergers). In general, the effect of MBH-MBH mergers is 
to decrease the spins of the most massive MBHs in the 
case of coherent post-feedback phase, while it increases 
their spins in the case of chaotic post-feedback phase. 
Broadly speaking the trends are in agreement with those 
found by Fanidakis et al. (2010, their "chaotic accre- 
tion") despite the different modeling and approach. In 
fact, in our model MBHs all accrete coherently during 
merger-driven episodes (thus increasing their spins) , but 
the post-quasar evolution (decline of the accretion rate, 
MC accretion, MBH-MBH mergers) broadens the spin 
distribution, overall decreasing spins at later times. 

In Fig. HI we focus on active MBHs. MBHs accret- 
ing at high rates, /Edd > 0.1, have very large spins 
at all 2; > 2. These are for the most part MBHs in 
the "quasar" phase. At 1 < z < 2 more systems are 
caught in the decline phase, and here, especially at high 
MBH masses, it becomes crucial whether accretion oc- 
curs mostly chaotically or coherently after the quasar 
phase (cf. left and right columns). This mass (> IO^Mq) 
and reshift (0.5 < z < 1) range is the most suitable to 
probe how feedback affects the angular momentum of 
nuclear gas. For low-mass BHs the spin distribution is 
mostly insensitive to the chaotic and coherent models 
(green triangles in Figl9|), while there is a stronger im- 
pact on high-mass holes (orange squares in Fig|9|). There- 
fore the changes most affect the high-luminosity end of 
the luminosity function. Finally, at z < 0.5 differences 
tend to disappear, as MC accretion becomes the domi- 
nant feeding mechanism. The behavior of lower accretion 
rate systems is similar, although more and more systems 
are in the "decline" phase, and lower spins become more 
common if accretion occurs chaotically during this phase. 

The spin distribution, however, does not necessarily 
map the radiative efficiency distribution (Fig. [TOl also 
note that the mean efficiency differs, mathematically, 
from the efficiency corresponding to the mean spin), if 
a large population of sources have geometrically thick 
and optically thin accretion discs where the radiative ef- 
ficiency may be suppressed with respect to the mass-to- 
energy conversion efficiency. Especially at low luminosity 
and low redshift, many AGN are radiatively inefficient 
sources and the MBH spin is not relevant in determin- 
ing their radiative efficiency. At low luminosity and high 
redshift AGN are a mixture of low mass MBHs accreting 
at high rates in radiatively efficient fashion, and higher 
mass MBHs accreting at low rates, while at low redshift 
most low luminosity sources are genuinely inefficient ac- 
cretors. At high luminosity, instead, the signature of spin 
evolution with redshift is evident. Fig. |9] and Fig. [10] 
show that our model is in very good agreement with the 
evolution of radiative efficiency with redshift derived ob- 
servationally by Li et al. (2012, see their Fig. 7), and by 
the suggestion of Li et al. and Shankar et al. (2011) that 
radiative efficiency may increase with black hole mass. 

Finally, in Fig. [11] we show the distribution of MBH 
spins before and after a MBH-MBH coalescence. As ex- 
pected, if MBHs have large spins pre-coalescence, MBH 
mergers tend to spin down the systems, and viceversa 
if MBHs have low spins pre-coalescence, MBH mergers 
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Figure 8. Spins in galaxies of different morphologies as a function 
of MBH mass (mean and 1-cr dispersion). Filled circles: gas-poor 
galaxies (spheroids). Stars: gas-rich galaxies (discs). Bottom: z < 
0.5. Top: 0.5 < z < 1. Lower mass MBHs in gas-rich galaxies tend 
to spin less rapidly than higher mass ones, and also less rapidly 
than MBHs in gas-poor galaxies. 

tend to increase spins. We here show the spin of the pri- 
mary MBHs in a binary prior to coalescence, as this is 
the quantity that gravitational wave observatories such 
as eLISA can measure. 

7. DISCUSSION AND CONCLUSIONS 

We have developed a model for the evolution of MBHs 
that takes into account several physical mechanisms of 
MBH growth: MBH-MBH mergers, merger-driven ac- 
cretion, stochastic accretion, and accretion of recycled 
gas. This model, however, does not include MBH feed- 
ing through disc instabilities, nor the burst phase of re- 
cycled gas feeding in elliptical galaxies. Under a series 
of plausible assumptions we have derived the growth of 
MBHs, the properties of the AGN population and the 
evolution of MBH spins. Our approach produces a pop- 
ulation of MBHs and AGN consistent with the observed 
one in terms of, e.g., luminosity function of AGN, re- 
lationship between MBHs and their hosts, high-redshift 
quasars. 

The main results of our models of MBH evolution can 
be summarized as follows: 

• At high-redshift MBHs grow mostly by merger-driven 
accretion, while at later times other channels become 
more important. In gas-rich galaxies, MC accretion dom- 
inates the growth of low-mass (< lO^Mg) MBHs at 
z < 2. In gas-poor galaxies, MBH-MBH mergers are 
the main growth channel, especially at high MBH mass 
(> lO'Me); 

• The mass of most active black holes decreases with in- 
creasing cosmic time, in "anti-hierarchical" fashion. Sus- 
tained accretion grows the most massive black holes since 
early times without overproducing the MBH population 
as a whole; 
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Figure 9. Spins in AGN, selected by mass and Eddington ratio. 
Triangles: MBH mass W^Uq < M < IO^Mq. Squares: MBH 
mass > IO^Mq. Top: QSOs accreting at high accretion rate (in 
Eddington units) f^dd > 0.1. Bottom: AGN accreting at all accre- 
tion rates f^dd > 0.001. Left: chaotic accretion during the decline 
of the quasar phase. Right: coherent accretion during the decline 
of the quasar phase. The signature of feedback on the angular mo- 
mentum of nuclear gas is strongest on high mass MBHs (compare 
right and left panels). In general AGN spins tend to decrease at 
late cosmic times. 
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Figure 10. Top: radiative efficiency of AGN of different luminosi- 
ties. Diamonds: 10*'' erg s~^ < Lbol < 10** erg s~^. Asterisks: 
^bol > 10** erg s-l. Bottom: spins of the same AGN. 30-40% of 
low luminosity sources (Lbol ~ 10*" erg s~^) are genuinely inef- 
ficient accretors. Sources with Lbol > 10** erg s~^ are powered 
by efficient accretors, and spin solely determines their radiative 
efficiency. 



Figure 11. Spins of the primary MBHs in a binary prior to co- 
alescence (pink filled histogram) and spin of the newly merged 
MBH post coalescence (violet hatched histogram) . Left: chaotic 
accretion during the decline of the quasar phase. Right: coherent 
accretion during the decline of the quasar phase. 



• MBH spins tend to be larger at redshifts z > 2, typ- 
ically a > 0.8. They result from massive, coherent ac- 
cretion events triggered by major mergers. This result is 
in general agreement with the trend found by Barausse 
(2012) using a complementary, more refined approach to 
modeling galaxy evolution; 

• A significant drop in the average value of MBH spins 
takes place at ^ < 2. This is caused by the increas- 
ing number of dry MBH-MBH mergers at lower red- 
shifts in the case of spheroids. Additionally, a dramatic 
drop is predicted at z < 0.5, for low-mass MBHs in gas- 
rich galaxies. This is due to low-mass, chaotic accretion 
events involving capture of molecular clouds; 

• In general, in gas-rich galaxies at 2 < 1 low-mass MBHs 
tend to spin slightly less rapidly than high-mass MBHs. 
The difference is less pronounced in gas-poor galaxies. 
As MBH mass increases the distributions become more 
similar. At the highest masses (> IO^Mq) the statistics 
are poor in the case of disc galaxies (between a couple 
and ~50 objects per bin), as the most massive among 
MBHs tend to reside in gas-poor galaxies. 

• If outflows do not affect the angular momentum of nu- 
clear discs, and accretion proceeds coherently both in the 
quasar and decline phase, the spin distribution and its 
evolution is not very differed for highly accreting vs low 
accreting MBHs. Differences are clearer in the popula- 
tion of the most massive black holes (> lO^Mg). 

• If quasar feedback disrupts the nuclear disc feeding the 
MBH and accretion proceeds chaotically in the decline 
phase, the spin distribution and its evolution shows a 
stronger dependence on mass, but not on morphology, 
again, except at the highest masses (> lO^Mg) and z < 
1. The same comment on statistical significance as above 
applies here as well. 

Qualitatively similar results have been obtained by Li, 
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Volonteri et al. 



Wang & Ho (2012) on observational grounds. They in- 
ferred the MBH spin evolution by tracing the evolution of 
the radiative efficiency of accretion flows, using the con- 
tinuity equation for the MBH number density. Both our 
and Li et al. results seem to contradict the predictions 
of the 'spin paradigm' scenario according to which the 
jet production efficiency - and therefore, radio-loudness 
of AGN - should reflect the MBH spin distribution and 
its evolution (Wilson & Colbert 1995; Hughes & Bland- 
ford 2003). Applying such a scenario to our results, 
one should expect the radio-loud fraction of AGN to be 
much larger at high redshifts than in the present epoch. 
At least in the case of quasars, an opposite trend has 
been inferred, i.e. such a fraction has been suggested 
to decrease with redshift (Jiang et al. 2007), although 
Volonteri et al. (2011c) flnd that the radio-loud fraction 
is roughly constant with redshift for the most luminous 
sources {L > 10^^ erg/s). For high- redshift blazars pow- 
ered by M > lO^M,;, MBHs, also, ac tivity seems to peak 
around z 4 iVolonteri et al.l (|2011l ) . Based on these in- 
dications, it may be that spin plays a more important 
role in powering jets at high accretion rates. Intermit- 
tent jet production in high accretion-rate AGN (Sikora 
et al. 2007), may occur in similarity to that observed in 
the galactic micro-quasar GRS 1915-1-105 (Fender et al. 
2004). However this does not explain the finding that 
radio-loud quasars reside in more massive and denser 
environments than radio-quiet ones (Shen et al. 2009; 
Donoso et al. 2010). Regarding low-accretion rate AGN, 
given the similar MBH spin distributions in the 'disk' 
and 'spheroid' sub-populations no double upper-bound 
pattern is expected to emerge in the radio-loudness vs. 
Eddington-ratio plots, in contrast to that found by Sikora 
et al. (2007, but see Broderick and Fender 2011 for a 
discussion on how large-scale environmental effects may 
affect low-frequency, low resolution radio loudness mea- 
surements, causing a larger dispersion in observed radio 
loudness) . 

It should be noted, however, that all these contradic- 
tions do not necessary jeopardize idea of the dominant 
role of rotating BHs in powering of extragalactic jets. 
This is because the Blandford-Znajek mechanism respon- 
sible for the extraction of rotational BH energy depends 
not only on the BH spin, but also on the magnetic flux 
threading the MBH. As Lubow et al. (1994) demon- 
strated it is rather impossible to cumulate in the center a 
sufficiently large magnetic flux to power strong jets if the 
accretion proceeds via geometrically thin disks. Hence 
the AGN that become efficient jet producers are those 
that have passed through the super-Eddington, advec- 
tion dominated accretion phase at least at the beginning 
of the massive accretion event, or those which were oper- 
ating in the very sub-Eddington, advection dominated 
accretion regime prior to the massive accretion event 
(Igumenshchev 2008; Cao 2011; McKinney et al. 2012). 
Therefore noting that radiogalaxies as well as radio-loud 
quasars reside on average in denser environments than 
radio-quiet AGN (Wake et al. 2008; Mandelbaum et 
al. 2009; Lietzen et al. 2011), this might imply that 
the broad range of AGN radio-loudness (see Sikora et 
al. 2007 and refs. therein) is mainly determined by the 
range of environmentally conditioned net magnetic fluxes 
collected in AGN centers. 
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APPENDIX 

BLACK HOLE GROWTH 

We here recall how one can describe the growth of a MBH as a function of a constant or variable Eddington rate. 
We start by defining f^dd = L/L-Edd, and ^Edd = Mc^/tEdd, where t^dd ~ i^cm ~ ^-^^ ^Edd represents 

the Eddington fraction. Therefore, if the accretion rate is Mi„, and M is the mass that goes into increasing the MBH 

L = eM„jC^ = /Edd^EddC^ (Al) 



mass: 



and dM = (1 - £)dMin and 



one obtains 



dM l-e 2 '^-<^r Mc^ 

— TT = ^EddjEddC — jEdd- , (A2) 

at € e lEdd 



dM _l-e 

If e and f^dd are constant over the time of integration, then: 



— = ^hddt^l^dt. (A3) 



J M-^dM^ j ^-If^^^t^l^dt, (A4) 

and the MBH mass grows as 

M{t + A<) = M{t) exp (^/Edd^i^^) , (A5) 
while if, for instance /sdd is a function of time, one has to self-consistently integrate: 

M-'dM^ f ^t^lJj,^^{t)dt, (A6) 



as shown in section 4.2. 



